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We present tight binding molecular dynamics simulations of Ceo collisions 



and with cells containing 1140 atoms. The results of our simulations are in 



^ ■ on the reconstructed diamond(lll) surface, carried out with an 0(N) method 
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f~^ ■ very good agreement with experiments performed under the same impact con- 

0^ . ditions. Furthermore our calculations provide a detailed characterization of 



the microscopic processes occuring during the collision, and allow the identifi- 



ed 
i-rt ■ cation of three impact regimes, as a function of the fullerene incident energy. 
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O ■ Finally, the study of the reactivity between the cluster and the surface gives 



insight into the deposition mechanisms of Cgo on semiconducting substrates. 
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Since the discovery [|l| of the unique stabihty of Ceo molecules against surface induced 
fragmentation, the investigation of fuUerene interactions with solid substrates has become 
an exciting field of research^"^. Collisions of Ceo on several surfaces have been reported 
to induce no fragmentation of the impinging molecules, for initial kinetic energies (Ek) up 
to a few hundred eV [|I|]. This property is unprecedented in molecular-ion surface-induced 
phenomena and has earned Ceo the name of resilient molecule. 

Of particular interest is the case of Ceo impacts on semiconducting surfaces which can 
form covalent bonds with the cluster. The study of these bonds is a challenging problem, 
from a fundamental point of view. Furthermore such an investigation is an essential pre- 
requisite for the understanding of fullerene deposition on non-metallic substrates. These 
deposition processes have recently been proposed as a valuable way of synthesizing thin 
films which would retain specific characteristics of the incident clusters. 

In order to understand the physics of Ceo collisions on solid substrates, it is crucial to 
reveal the processes occurring at the surface. Whereas these processes are difficult to probe 
experimentally |I]], they can be investigated by computer simulations. Recently, simulations 
of Ceo impacts on surfaces (hydrogenated diamond and graphite) have been performed by 
using empirical classical potentials 0]. These calculations constitute a useful starting point 
for the study of fuUerenes on surfaces. However, a more accurate description of the cluster- 
surface interactions is desirable. In principle this can be obtained by molecular dynamics 
(MD) simulations in which the interatomic forces are derived at each step from quantum 
mechanical (QM) computations. In practice, realistic simulations of Ceo interactions with a 
surface require cells with at least several hundreds of atoms. Until recently QM simulations 
of systems of this size have not been possible because of the computer time required by 
conventional algorithms, which grows as the cube of the number of atoms in the system 
(N). Indeed the only semiempirical QM study of Ceo collisions with a surface (hydrogenated 
diamond) reported to date 0] has been limited to small MD cells. 

This situation has changed with the introduction of methods [|Ty-[T2| for electronic struc- 
ture calculations and MD simulations, which are based on algorithms whose computational 



workload grows linearly with the system size. These approaches (often referred to as 0(N) 
methods) have opened the way to QM-MD simulations of systems much bigger than pre- 
viously accessible, and therefore to the study of problems usually out of reach for QM 
calculations. 

In this Letter we present QM-MD simulations of Cgo impacts on a semiconducting surface 



as a function of the molecule incident energy, carried out with an 0(N) method |[TT|, and 



with MD cells containing 1140 atoms. Interatomic forces are described according to a tight 



binding (TB) Hamiltonian [^, derived from first principles calculations. We considered the 
clean (2 x 1) reconstructed diamond (111) surface. This surface can be expected to form 
bonds with the impinging Ceo, since its uppermost layer contains three fold coordinated 
IT bonded atoms; it therefore represents a good candidate for the study of covalent bonds 
between the fuUerene and a semiconducting substrate. Our work constitutes the first QM 
study of such interactions and the first large scale TB-MD simulation based on an 0(N) 
methodology. The results of our calculations are in very good agreement with experiments 
carried out under the same impact conditions [I|-|^. Furthermore they allow us to charac- 
terize in detail the microscopic processes occuring during the collision and thus to identify 
three different impact regimes as a function of Cqo incident energy. 

In our calculations, we used the 0(N) method of Ref. [|ll|. This approach is based on 



an energy functional with implicit orthogonalization constraints and on a localized orbital 



(LO) formulation. We adopted the TB Hamiltonian {"H) proposed by Xu et al. [[13| , p!4| . The 
LOs were centered on atomic sites, extending up to second neighbors (in the notation of 
Ref. [jll|, Nh = 2) [|15|. We simulated neutral Ceo molecules colliding with a reconstructed 



C(lll) surface at initial energies Ek= 60, 80, 120, 150, 180, 210, 240, 300, 400 eV. Each 
of the nine runs lasted from 0.3 to 1 ps. The efficiency of the 0(N) algorithm allowed us 
to carry out all calculations on workstations. Cgo impinged upon the surface at normal 
incidence, oriented in such a way that a double bond between two hexagons faced a Pandey 
chain of the surface. The C(lll) substrate was represented by a slab composed of twelve 
layers, each containing ninety C atoms. The dimensions of a layer were 22.72 x 21.86 A in 



the X and y direction, respectively. The slab was terminated on each side by a reconstructed 
surface, and periodic boundary conditions were applied along x and y. Nine surface layers 
were allowed to move during the simulations. According to the size of our MD cell and to 
the number of mobile layers, we estimate that the shock wave produced by a collision of Ceo 
with the substrate is echoed back to the impact region about 0.2 ps after the impact. This 
is the time interval during which we observe the formation of bonds between the cluster 
and the surface. Therefore the characterization of the three impact energy regimes (ERs) 
proposed in the following is not affected by finite size effects. 

Before starting the simulation of Ceo collisions on the surface, we optimized the slab and 
the molecule geometries independently. Our results for the surface reconstruction compare 
well both with those of a conventional (extended orbital [EO]) calculation using the same Ti 



17 1, and with the findings of first principles computations |T8|. For example, we find that 



the average bond distance in the surface chains is 1.44 A and that the change in bond length 



between the second and third layer is 6%; the corresponding values reported in Ref. |T^ 



are 1.48 A and 6%, whereas those of Ref. |T^ are 1.44 A and 8%. However our calculation 



seems to overestimate the dimerization of the surface chains (4% versus 0.2% in Ref. |T7 



and 1.4% in Ref. [T^). The two characteristic bond lengths of Cgo calculated within the 
LO formulation differed by 5% (single bond) and 2.5% (double bond) from those optimized 
with EO [l^ . Also in this case the TB results [^] are in good agreement with those of first 



principles calculations [^ 



We now turn to the discussion of the different regimes of Ceo-surface interactions which 
we observed in our simulations. For Ek=60, 80 and 120 eV (low ER) the molecule does not 



form any bond with the surface |jT6|. The minimum distance between Ceo atoms and the 
substrate (D) is approximately constant during the collision and very close to 2 A, as seen 
in Fig. 1, where we show D as a function of the simulation time (t). The impact provokes 
large distortions in the fuUerene cage: the cluster height is decreased from 7.0 to 3.9, 3.4 
and 3.lA , for Ek= 60, 80 and 120 eV, respectively, and the flattened molecule adapts its 
shape to that of the reconstructed surface (see Fig. 2). For Ek=120 eV, a few bonds of the 



molecule are broken during the collision, whereas no bond breaking is observed at lower 
energies. This can be seen in Fig. 3 which displays the number of differently coordinated 
sites of Ceo as a function of t. In all of the three cases, after the impact the cluster leaves 
the surface with its original shape, without any defect. The kinetic energies of the centre 
of mass (Ecm) of the outgoing molecule is shown in Fig. 4. In this ER, E^m shows a clear 
dependence upon Ek. In particular Ecm is proportional to Ek for incident energies up to 80 
eV, i.e. in the energy interval in which no bonds of the cluster are broken. 

For Ek larger than 120 and smaller than 240 eV (medium ER), the molecule does not only 
suffer large distortions upon impact on the surface, but also forms bonds with the substrate. 
During the collision, some bonds of the Cqo cage are broken, resulting in the formation of 
twofold coordinated sites. The presence of defects in the Ceo structure causes the molecule 
to react with the surface, at variance with what happens in the low ER. Indeed, after some 
oscillations around 2 A, D suddenly decreases to a value of about 1.5 A, i.e. to a distance 
typical of sp^ like bonds (see Fig. 1). Not surprisingly the surface atoms involved in the 
bonding belong to the topmost layer, containing three-fold coordinated sites in the Pandey 
reconstruction. In this ER, the number of bonds between Ceo and the surface depends 
weakly on Ek and, during the collision, is approximately two. For Ek=150 eV, we find that 
these bonds are stable after the impact and that the cluster is adsorbed on the surface (see 
Fig. 2). We note that a decrease in scattering intensities around 140 eV has been found in 
collision experiments on a graphite substrate p|,0, and that adsorption of intact Cgo was 
invoked by Busman et al. p[ in order to explain this observation. At Ek=180 and 210 eV, 
the Ceo-surface bonds break after the collision and the molecule leaves the substrate with 
Ecm = 8 and 15 eV, respectively. Our findings compare well with those of Ref. 0, where the 
Ecm of fuUerene ions scattered from graphite is found to be about 10 to 20 eV and nearly 
independent of the impact energy, for Ek larger than 140 eV. In this impact regime, after 
the collision, the topology of the molecule is still that of a cage, although defects are present 
with respect to the original shape (see Figs. 2 and 3). 

Upon impact at Ek > 240 eV (high ER), we observe the formation of several bonds 



between Ceo and the surface. Their average number increases with the incident energy, 
varying from about six to sixteen. Shortly after the colhsion, many bonds break within the 
molecule and the cage structure of the fuUerene can no longer be identified (see Fig. 4); a 
disordered structure is formed and adsorbed on the substrate. After the 300 and 400 eV 
collisions, some of the Ceo atoms form twofold coordinated chain structures (see Fig. 2); the 
formation of such chains was observed also in MD studies of fuUerene melting |2^ . 



Whenever Ceo left the surface, we did not observe any fragmentation of the molecule. 
However a rupture may be expected to occur some time after a collision if the internal energy 
(Ein) of the bouncing molecule (see Fig. 4) exceeds its stability threshold. MD investigations 



of fuUerene melting |2^ have shown that Ceo becomes unstable when its internal kinetic 
energy (Ejn) is between 30 and 40 eV. 

In the low and medium ER the topology of the reconstructed surface is unaffected by 
the collision, although a considerable fraction of Ein is transferred from the molecule to 
the surface (e.g. as much as 50% for E^ = 180 and 210 eV). In the medium ER only the 
three fold coordinated atoms of the topmost layer bond to Ceo atoms (see Fig. 2). In the 
high ER the fourfold coordinated surface atoms (belonging to the second layer) are also 
involved in bonding to the molecule. However they remain fourfold coordinated since their 
bonds with the third layer atoms are broken. This bond breaking induces a local surface 
dereconstruction from a tt bonded chain geometry towards a (1 x 1) ideal arrangement (see 
Fig. 2). 

We finally note that the connection between Ceo enhanced reactivity with the surface 
and the presence of defects in the molecule might characterize deposition processes of the 
fuUerene on other semiconducting surfaces. Therefore the variety of STM images produced 
by covalently bonded Ceo on, e.g., Si(lll) 0, might correspond to molecules with different 
types of defects. 

In summary, we have identified three impact energy regimes of Ceo collisions on C(lll). 
For Ein < 120 eV, the molecule bounces off the surface without ever forming bonds, and 
recovers its original shape after severe distortions. In a second regime. Ceo forms covalent 
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bonds with the surface. The cluster can be either adsorbed on or leave the substrate, 
with defects in its original structure. The energy at which we observe absorption (150 eV) 
corresponds to the energy at which a decrease of fullerene ions scattered from graphite is 
observed experimentally. In the third regime (Ek > 240 eV ) the molecule breaks after the 
impact and pieces of the broken cage form stable bonds with the substrate, which can induce 
a local surface dereconstruction. The formation of bonds between C(lll) and the fullerene 
is always accompanied by the formation of defects in the molecule. When the molecule 
bounces off the surface, the kinetic energy of its centre of mass shows a clear dependence 
on its incident energy in the low ER; on the contrary Ecm depends weakly on Ein in the 
medium ER and acquires values in agreement with those measured under the same impact 
conditions. Finally we have shown that QM simulations for systems containing thousands 
of atoms, so far not possible with conventional methods, are now feasible by using 0(N) 
approaches. In the simulations reported in this work, the gain in computer time with respect 
to conventional O(N^) methods is estimated to be of the order of 1000. 
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FIGURES 
FIG. 1. Minimum distance (D) between Cqo and the surface as a function of the simulation 

time: soUd, long dashed, dashed and dotted lines correspond to Ek= 120, 150, 180 and 210 eV, 

respectively. 

FIG. 2. Snapshots of the uppermost five layers of the slab and of Cgo for different E^ and at 
given simulation times (t): Ek = 120 eV and t=0.14 ps (upper panel); E^ = 150 eV and t=1.00 
ps (middle panel); E^ = 300 eV and t=0.33 ps (lower panel). (See text). Red and green spheres 
denote threefold and fourfold coordinated atoms of the slab and black spheres indicate Cqo atoms. 

FIG. 3. Twofold and fourfold coordinated Cgo atoms (Nat) as a function of the simulation 
time for selected E^: dotted, solid and dashed lines correspond to E^ = 120, 180 and 300 eV, 
respectively. 

FIG. 4. Final energies (Ef ) of Cgo bouncing on the surface as a function of E^: diamonds and 
crosses indicate the kinetic energy of the centre of mass and the internal energy of the molecule, 
respectively. 
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